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ABSTRACT 

We show that the Husain-Kuchaf model can be described in the framework of BF 
theories. This is a first step towards its quantization by standard perturbative QFT 
techniques or the spin-foam formahsm introduced in the space-time description of 
General Relativity and other diff-invariant theories. The actions that we will consider 
are similar to the ones describing the BF- Yang-Mills model and some mass generating 
mechanisms for gauge fields. We will also discuss the role of diffeomorphisms in the 
new formulations that we propose. 
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I Introduction 



BF theories were introduced by Horowitz |1|] and, independently, by Blau and Thomp- 
son in the late eighties and early nineties. They are simple diff-invariant models 
that have been extensively studied as a testbed for quantization techniques in the 
absence of metric backgrounds. They have also been used as building blocks for 
physical theories such as gravity and Yang-Mills [^. From the point of view of 
their quantum treatment it is interesting to point out that they can be quantized both 
by using traditional Quantum Field Theory (QFT) techniques and methods specially 
tailored to deal with diff-invariant background-free actions such as spin networks and 
spin-foams 0. 

The Husain-Kuchaf model (HK) |^ is an interesting diff-invariant system with 
local degrees of freedom that actually mimics many of the features of 3+1 dimensional 
general relativity (GR). In fact, the solutions to the Einstein equations (with an 
internal 5*0(3) group) can be seen as a subset of the solutions to HK. Its Hamiltonian 
description is very similar to the Ashtekar formulation for GR because both share the 
same phase space and most of the constraints -only the Hamiltonian constraint is 
missing in HK-. Its absence means that the model describes equivalence classes of 
metrics^ in the spatial slices of a 3+1 foliation of space-time, but no time evolution. 
It is widely believed that a successful quantization of HK, and the implementation of 
efficient computational techniques to deal with it, would be important steps towards 
the quantization of full GR. 

The standard action for the HK model is very similar to the self-dual one for 
GR because they essentially differ only in the internal gauge group (5*0(3,1) for GR 
and 5*0(3) for HK^). This explains the similarities of their Hamiltonian descriptions. 
Some other action principles for this model have also appeared in the literature M- 



10|. They provide different Hamiltonian formulations of the theory and solve some 



apparent shortcomings of the original formulation -such as the description of non- 



•^Under both difFeomorphisms and SO{3) internal rotations. 

''Actually the self-dual action in the Euclidean case can be obtained from the HK action by 
adding a term to it as shown in [BI . 
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degenerate 4-metrics and the possibility to couple ordinary matter-. One of them 
(the two-connection formulation Q]) will be used in this paper as the starting point 
to obtain the HK model by coupling BF Lagrangians. The action that we give is very 
similar in form to the ones appearing in |^ for BF- Yang-Mills (BFYM) models, in 
TT| for a mass generating mechanism for gauge fields, or in P to discuss the problem 



of infrared divergencies in the spin-foam quantization of the BF theory. 

There are several reasons that lead us to believe that the actions that we consider 
here are of interest. First of all they are structurally very simple because they consist 
only of several BF terms plus quadratic interactions. At variance with the known 
formulations, they have quadratic terms that may make it simpler to use perturbative 
QFT techniques. Also, we want to emphasize that we are not resorting to the usual 
(and rather trivial) trick of taking a BF model and impose restrictions on the 2- 
form field B with Lagrange multipliers because we want to preserve the quadratic 
character of the action as much as possible. Second, BF theories have also been 
quantized by using the so called spin-foams that provide a picture of quantum space- 
time geometry. As is shown in one may need to add a "cosmological constant" 
term to remove apparent infrared divergencies that pop up in the computation of 
transition amplitudes^. This term is very similar to the ones that we introduce here 
so we hope that the strong similarity of our formulation with this one will provide a 
completely different way to deal with the HK model. Some of the issues that are not 
yet fully understood in the spin-foam framework, such as the relationship between 
the presence of local degrees of freedom and triangularization independence may be 
illuminated by using the actions that we give here. 

An additional reason to consider different formulations for old models is the fol- 
lowing. As it has been discussed by some authors [|12|, [jl3[ the non-renormalizability 
of GR can be traced back to the non-invertibility of the quadratic part of the ac- 
tion after gauge fixing of its symmetries. In the case of the Hilbert-Palatini action 
there is, in fact, no quadratic term to build field propagators. Something similar 



^The functional integrals in Q are performed with the naive measure in field space without 
introducing the terms needed to take care of the second class constraints present in the theory. 
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happens for the previously known actions for the HK model because all of them lack 
a quadratic part. A possible way out of this is to consider the non-quadratic actions 
as interaction terms of theories modified by the addition of quadratic kinetic terms 
and try to extract some information from these modified models. For this approach 
to work it must be possible to add these kinetic terms in a consistent way, i.e. the 
fully interacting theory must be a consistent deformation (in the sense of |T^) of its 
quadratic part. As it turns out it is actually impossible to modify the theory by 
adding kinetic terms (both diff-invariant and in the presence of a metric background) 
to have a consistent formulation 1113 where the Hilbert-Palatini or the HK actions 



appear as interactions. This is a consequence of the fact that these actions are built 
in terms of 1-form fields, so it is interesting to study whether new formulations in 
terms of different objects -such as 2-forms- will allow us to successfully incorporate 
consistent kinetic terms, or hopefully, define propagators. 

A final issue that we will discuss is the role of diffeomorphisms in HK actions. 
We will show that, as in the 2+1 gravity case (a BF theory itself), one can actually 
write down non-diff invariant actions leading to the same field equations and gauge 
symmetries (on shell) as those of the standard formulations. 

The paper is organized as follows. After this introduction we give, in section 
H, the BF description of the HK model -following a brief discussion of some of the 
other known formulations- and the general mechanism at work here and in other BF 
descriptions of known theories. Section HI contains our conclusions and comments. 
We end the paper with an appendix where we discuss the detailed Dirac analysis for 
the actions that we present in the paper. 

II BF Formulation of the Husain-Kuchaf Model 

BF theories in four dimensions are described by the action 

Sbf= I d^xfi'^^'^'^B^uFi,, (1) 

J M 

where fj"'^'^'^ is the four-dimensional Levi-Civita tensor density, F^^, is the curvature 
of a gauge connection 1-form (taking values in the Lie algebra of a certain gauge 



group G) defined as = 2d[aAl-^ + [Aa, Ab]\ and Babi is a 2- form defined in the dual 
Lie algebra. Here and in the following we use tangent space indices a, b, c,. . . In this 
paper we will restrict ourselves to G = 5*0(3); and use internal indices i, j, k, ... = 
1,2,3. The totally antisymmetric Levi-Civita tensor will be denoted as eijk- The 
differentiable manifold A4 is taken to have the topology = IR x S with E a 
three-dimensional compact manifold without boundary. 
The action for the Husain-Kuchaf model is 

Shk = I d'x r''%ke\elF^',. (2) 

where F*^ = 2d[aAl-^ + e^^^ AajA^k is the curvature of a S'0(3) connection A\ and e\ are 
three 1-forms (notice that they are not proper tetrads because there are only three of 
them). As shown in it can be rewritten (modulo surface terms) as 

S'hk= I d'xfi^'^'huhd. (3) 

JM 

12 1.1. 

where F au and F au are the curvatures of two 5*0(3) connections = 29[a74^] + 

...1 1 2. 2. -w^ 2 

e"-^ Aaj Abk and F^^ = 2d[aAl-^ + e"-^ Aaj Abk- This can be easily seen by defining 

2 1 

Cai = Aai — Aai and using the fact that 

Kb = Kb + 2V[ae^] + e'^'^ea^Cbk (4) 

(the covariant derivative acts on internal 5*0(3) indices as Va-^i = daXi+e^-' AajXk 
and is defined in an analogous way). By using (|) we can write the action (^ as 

Shk= I d^xf]'''''''\hb^Kd + '^KbWced^ + e''''ea^eb,hdk]] (5) 

dropping the first term due to its topological character and the second after integrat- 
ing by parts and using the Bianchi identities we end up with the action (H). Notice 

1 . 2 . 

that both F\^ and are defined in the Lie algebra of 5*0(3) so we must use the 
invariant metric 5"^^ to build the action. The field equations coming from (|^) have the 
beautifully symmetric form 

(6) 



Let us consider now the following action^ 

r rll 22 12n 

Sbfhk = I d X fj" " BabiFl^+ BabiFl^+ BabiB'^d ) (7) 
JM '- -' 

1 2 . 

obtained by coupling two BF theories via a rj"'^'^'^ B^bi B^^ interaction term. The field 
equations coming from are now 



Fl, + Bl, = r V[a 5L] = 

2 . 2 (8) 



The equivalence between these equations and (H) is obvious. The new action Sj^p^^^ 
is both S'0(3) and diff- invariant. To understand how we get local degrees of freedom 
from the two topological BF models it suffices to realize that an action obtained by 
adding two different BF terms has two independent sets of symmetries that are par- 
tially broken by the rj Bau B\^ term; specifically, without this coupling the action 
would be invariant under two independent sets of 5*0(3) rotations and diffeomor- 
phisms. The reduced gauge symmetry accounts for the new local degrees of freedom 
described by (0). By using the terminology introduced by Smolin in the action 
(0) describes the HK model as a constrained topological field theory. 

There is an evident, one to one, correspondence between the solution spaces to 
equations (H) and (|) so, classically, both actions (|]) and (|^) describe the same theory. 
The equivalence of their quantum formulations is, however, less obvious. As Sbfhk is 
quadratic in i?^^ and B\^ one would be tempted to say that the effect of performing 
the functional integration in these fields would be equivalent to substituting the field 
equations for them in the action, in which case we would end up with the 2-connection 
action (^) and a functional integral measure involving only the connections. However 
this is not correct because the functional measure is not the naive one but must be 



modified 16 to take into account the existence of second class constraints in the 



Hamiltonian formulation derived from (|^). An intuitively clear way to understand 

^Here we are following a suggestion due to Giorgio Immirzi to couple two Yang-Mills theories by 
using their BFYM formulations. 



this fact is to realize that, as a consequence of the Bianchi identities, only those 

1 . 2 . 

components of -B^^ and that cannot be written as covariant derivatives (in terms 

1 . 2 . 1 . 2 . 

of the connections AJ^^ and A^^) couple to the curvatures and respectively. 

One must, in fact, remove them form the path integral in order to avoid problems 

such as the violation of the Ward identities that one encounters in BFYM models 

1 . 

m. After modifying the measure to circumvent this difficulty the integration in 

2 . 

and -B^j, is no longer trivial. The Dirac analysis of is rather involved because one 
expects to find many second class constraints. Most of them appear as consistency 
conditions for the solvability of the Lagrange multiplier equations that one must 
introduce in the total Hamiltonian. A way to partially alleviate this problem is to 
use Stiickelberg's procedure in by introducing auxiliary fields 77^ and in 
the action and trading some second class constraints for first class constraints. So we 
will also consider 



M 



(9) 



The field equations derived from (01) are 



[ F\, + 5^, - V[a^ 



b] 



2 k 





0, 



(10) 



0. 



From the first pair we find 



r Bi, = - Fl, + ^,,7^ 



B. 



ab 



ab 



"WlaV 



b] ' 



(11) 



so that the remaining equations give (|^). After solving them we can obtain and 

2 . 1 • 2 • 

-B^^ (modulo the arbitrariness in r/^ and from (|11]). Notice that the third pair of 
equations in (|T0|) is just a consequence of the first and the Bianchi identities. 

In addition to the S0{3) and the diff-invariance of (J^) this action is invariant 
under the transformations 



r SAi = 



6Ai = 



r 5Aj. = o 



(12) 



where \ and \ are arbitrary gauge parameters. A possible way to ensure that 

these are the only symmetries of the proposed action (and a first step towards the 

quantization of the model) is to study its Hamiltonian formulation. This is done in 

the appendix. One can see that the model is, in fact, equivalent to the 2-connection 

formulation after performing a suitable gauge fixing of the symmetries introduced by 

the Stiickelberg procedure. Here we follow a different path to prove the equivalence 

of Sbfhk and 5*^^^^^ with (||) at the quantum level by using the covariant symplectic 

techniques introduced in 

To this end we must look at the symplectic form in the space of solutions to the 

field equations This is obtained in a two-step process; first we write a closed 

(but degenerate at this stage) 2-form in the space of fields. For a Lagrangian of the 

type L(ti7, dw) depending on a s-form field w and its exterior derivative it is given by 

r dT 

Q= I dlu7Adl-r— , (13) 



dtu Ad- , 
s odzu 



where we define ^r- according to 



dL 



Lidw + de) = Lidw) + de /\ — h higher order . 

odzu 



Here we use the differential form notation to make a clear distinction between space- 
time and phase space objects. We must distinguish between the ordinary exterior 
differential in the spacetime manifold Ai (d) and the exterior differential in the field 
space(cll). In the same way we must make a distinction between the wedge product 
in both cases (A and A respectively). In the following we will avoid this by explicitly 
writing tangent space indices. Notice that d is defined in Ai but the integral in (0) 
is three-dimensional. 

We must now pull-back Q to the space of solutions to the field equations. When- 
ever this can be explicitly done (in those few cases where the solutions to the field 
equations can be completely parametrized [|T^) one gets a closed and non-degenerate 
2-form in the space of solutions to the field equations. The non-degeneracy implies 
that no gauge freedom remain^ and the objects that appear in the symplectic form 
parametrize physical degrees of freedom. For the action (|^) we have^ 

n = 2 [ d^xfi"''^ hAUx) AdBbJx) + (IAUx) AdBbdix)] , 

where 

~ahc -g 

the three dimensional Levi-Civita tensor density in the a;-coordinate 

1 . 2 . 1 . 2 . 

patch and A^, and B^^ are the pull-backs of the corresponding four dimen- 
sional objects onto S. After partially pulling back to the solution space by using the 

1 2 

field equations expressing B and B in terms of the curvatures we get^ 



^] = -4 



J^d'xfi''''' {dlil(x) h[^b - Wb]dA,i{x)} . (14) 



This is not yet the symplectic structure in the solution space because we still have 
to consider the remaining field equations but it coincides with the one derived from 
(|^). This not only proves the equivalence of the classical theories but also of their 
quantum formulations because their reduced phase spaces coincide and have the same 
symplectic structure. For the action (^ we obtain in a similar fashion 

n = J^d^x ff"" {2dl A\{x) Ad h^{x) + 2dl A\{x) Ad Bbci{x)- (15) 



^This is equivalent to the reduced phase space of the usual Hamiltonian formulation. 
^Notice that the A{x) and B{x) fields are time dependent; here d^x denotes the measure in E. 
^This coincides with the result in This example illustrates the efficiency of covariant sym- 
plectic methods as compared to the more traditional approach of Dirac. 



and, as before, after partially pulling back to the solution space by using the field 

1 2 

equations expressing B and B in terms of curvatures we get 

Vl = jj.^xfj"^^ {<lh\{x)h<lhci{x) +<l^\{x)h<lhci{x)+ (16) 

+2dlii(x)Adl['^,r/,,(a;) - hci{x)] +2dli^(x) Adl['^,f),,(x) - hci{x)]] 

which, after a little algebra gives precisely (|l^ and, hence, we also prove the classical 
and quantum equivalence of (|) with the actions considered before. 

It is interesting to compare this situation with some other similar models. Let 
us consider first the BFYM theory that is known to be amenable to perturbative 
treatment ||^. The BFYM action can be written as 

S,^,M. B] = Sym[A] - I d'x {Bl, + -F'J{B,a^ + -F.^dv^'ri'"' , 

jm g g 

where S'yAf[A] is the usual Yang-Mills action for A with a coupling constant g and 
7]"'^ is a background Minkowski metric. The term {B + F/ gY vanishes on shell thus 
making the BFYM and the YM models equivalent^. This term yields a propagator 
for B. Moreover Sym has, after a suitable gauge fixing, a well defined propagator for 
the Yang-Mills field A. 
In the BFHK case we can also write 

^2, B]= f dW'^'i Kb + Kb)ikd^ + Fed.) - S'hkIAi, A2] , (17) 

2112 

but now S'fjj^ does not have a free part and the quadratic form {B + F){B + F) is 
always degenerate in the A, B fields, even after gauge fixing. This prevents us from 
using standard QFT techniques because we do not have a suitable propagator even 
though we have a quadratic piece in the action. However (|l3) suggests a possible way 
out. The main problem that we face is related to the diff-invariance of the quadratic 

terms derived from Sbfhk- Looking at the mechanism that is working in (|T7|) , one can 

^''This argument, based on the resolution of algebraic field equations back in the action, relies on 
the following fact: Sbfym action generates a covariant phase space that is in one-one correspondence 
with the Yang-Mills one. Moreover, it induces a symplectic two form over the space of A, B fields 
that coincides, when restricted to the solution space, with the symplectic Yang-Mills form. 

9 



devise a straightforward generalization where the diff-invariance is broken through 
the use of a metric background g, so let us consider, for example, the action 

S" = S'nK- I + * Kb) (* k<i^ + hd^v'^'v'" ■ (18) 

JM 

where is the Hodge dual F^^. Once again, the equivalence with the HK model 
is due to the vanishing of the i?-term of the action in the covariant phase space. 
Clearly, although 5"' is not difF-invariant, the diff-invariance is recovered on shell. 
This formulation has the same perturbative objections that Sbphk- However, the 
presence of a metric background opens up the possibility of adding to the action S" 
kinetic terms for A and B. Then S" might be considered as an interaction term of an 
action with well-defined propagators. Of course this would be possible only if these 
(local) kinetic terms can be incorporated in a consistent way, both from the classical 
and quantum point of view. 

In addition to the BFYM models some other theories of this type have been 
studied in the literature. An interesting example is the mass generating mechanism 
for gauge fields proposed by Lahiri in |TT| and described by the action 



JM L 6 4 2 



(19) 



with Habc = 3V[a-B6c] and m is a mass parameter. The essential difference between 



our model and Lahiri's is related to how the B field appears in the action. In 
after disregarding total derivatives, one can write for the linear theory 



(20) 



but now the {dB — m * A^ term does not vanish on shell but it is proportional to 
the gradient of the gauge parameter for the A field. Another property of that 
makes his model closer to S", is the impossibility of using standard perturbative 
treatments |T9[. This is so because the free Lagrangian ( pO]) has two independent 



gauge symmetries 



6,Ai = daA' 62AI = 



10 



for arbitrary 0-forms A* and 1-forms while the full model has only a SO (3) symme- 
try. This problem is exactly the same that afflicts the Palatini gravitational theories 
in four space-time dimensions [jl5[. The open question is, again, the existence of 
consistent theories, in the perturbative sense, for interacting 1-form and 2-form fields 
implementing this type mechanism. As a systematic study of quadratic actions for 1- 
form and 2-form fields has not been carried out to date we hope that suitable theories 
can be found. 

Ill Conclusions and Comments 

As we have shown in the paper the Husain-Kuchar model can be obtained by coupling 
BF theories. The mechanism at work in the actions and (|^) is different from the 
standard procedure of writing a BF action and imposing constraints on the 2-form 
field with Lagrange multipliers. In the examples that we consider in this paper the 
local degrees of freedom of the model (three per space point) appear due to the 
breaking of some of the symmetries of the BF terms induced by the introduction of 

1 2 

a "cosmologial constant" term coupling the 2-form fields B and B. The main reason 
why we do this is that we want to have a well defined quadratic part in the action 
in order to define propagators (after a suitable gauge fixing of the symmetries) and 
treat the model by using standard perturbative QFT techniques. 

Of course, even in the presence of a quadratic part, the theory may not have suit- 
able propagators. This happens if the symmetries of the quadratic part do not match 
the symetries of the full action (i.e. the full action is not a consistent deformation of 
its quadratic part). In the present case we do not get a perturbatively quantizable 



action because, as shown in [|T^, its kinetic term is diff- invariant so it describes no 
local degrees of freedom and hence it has more symmetries than the full action. Any- 
way, it may be possible to modify these actions by the addition of local kinetic terms 
that may render the quadratic part invertible in such a way that the full action is a 
consistent deformation of it. If this can be achieved it may be possible to effectively 
recover the HK model in some limit, or at least extract useful information from its 
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behavior as an interaction term. One of the points that we emphasize in the paper 
is precisely the fact that writing the action in terms of new types of fields it may be 
possible to have kinetic terms not available in the known formulations. 

We have shown that one can actually derive the HK models from non diff-invariant 
actions. This is similar to GR because, as is well known ||2^, it is possible to derive 
the Einstein equations from a non-diff-invariant action, quadratic in the Christoffel 
connections, obtained by integrating by parts the second order terms in the expression 
for the scalar curvature. In this case the loss of symmetry is compensated by the 
appearance of many second class constraints that reduce the number of physical 
degrees of freedom to the two physical ones. Diff-invariance is kept only on-shell. 
Again, the availability of a background structure may permit the introduction of 
kinetic terms that could not be written without its help. 

We want to point out that covariant symplectic techniques are of great help to 
study the physical content of field theories without having to rely on the more tradi- 
tional Dirac analysis. We have shown that it is very easy to prove the classical and 
quantum equivalence of the actions that we have discussed here. This requires some 
comments. In principle, two non-trivially different actions]^ may lead to the very 
same field equations. This does not mean, however, that they are equivalent from the 
quantum point of view. To prove their quantum equivalence the actions must endow 
their solution spaces with the same symplectic structure. This does indeed happen 
in this shown above but this should not happen in general. 

In a precise sense it is not necessary to go through the painful procedure of using 
perturbative techniques as in to show that the BFYM action is equivalent, from 
the quantum point of view to the usual Yang-Mills one. Perturbation theory, if 
properly done (and disregarding difficult issues of convergence), reduces to canonical 
quantization. If one can argue, as we did before, that the solution spaces of the usual 
Yang-Mills formulation and BFYM coincide and are endowed by the action with the 

same sympletic structure their quantum formulations must coincide (which is the 

^^By non-trivial we mean that they not differ in total derivatives or cannot be transformed into 
one another by simple redefinitions of the objects in terms of which they are described. 
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result that one finds out, at the end of the day). 



IV Appendix: Dirac Analysis of the Coupled BF 
Action for the HK Model 

We give a detailed Dirac analysis for the action (|^) to show that we arrive exactly at 
the same conclusions that we got by using covariant symplectic techniques. We start 
by introducing a foliation of the space-time manifold = IR x E (where S denotes a 
compact 3-manifold without a boundary) defined by a scalar field t and a congruence 
of curves nowhere tangent to the t = constant hypersurfaces parametrized by t. In 
this way, if t is used as a "time coordinate" we have dt = Ct, i.e. time derivatives 
are given by the Lie derivatives along the tangent vectors t to the curves of the 
congruence parametrized by t. In the following we make no distinction between three 
dimensional and four dimensional tangent indices as it will be clear from the context 
if we are considering 3-dim or 4-dim objects. By using the identity fj"''"^^ = 4t[«^M^ 
where fj'^'^'^ is the Levi-Civita tensor density in three dimensions, and defining 



2 1 
~ahc TD i — rja 



V ^ bci = -T i 



1 1 



t- A\ ^ A\ 



t- A\ ^ A\ 




(21) 



we find 



(22) 



\h\Na[fl''''^bh'.-k) + 
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1 2 1 



The momenta canonically conjugate to the variable^ A^, A^, B^;, 5", A*, A*, 

9 19 
121' 2' 121-2-1212'-~ 

/^?, Vt^ Vt, Vl and r^* will be denoted as jr^, ^i, ^i, i>-, i>-, ^f, V'f. 

1 2 

and t/ij. They satisfy the following primary constraints 



1 



k = o 






1 

Tfi 

2 



2 






: 
: 








k 





0, 



(23) 



The Hamiltonian is 



+ 



+ 



(24) 



+ ^0 



+ 



and the total Hamiltonian 



7/" + n* 75"+ 



(25) 

1.212 ^ 2 1-2. 1.21 2 

where l\, I i, I i, m", Sj, and are Lagrange multipliers. By 

imposing stability of the primary constraints under the evolution defined by the total 
^^Functional derivatives with respect to Lt of the corresponding fields. 
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Hamiltonian Ht we get secondary constraints and equations involving the Lagrange 
multipliers. The secondary constraints are 



-0 

(26) 

:0. 



We also find a set of equations that the Lagrange multipliers must satisfy. In this 
case the stability of the secondary constraints does not give any new ones but only 
additional equations for the Lagrange multipliers. These equations must be treated 
with care. In many cases it is possible to show that they can be solved for any 
configuration of the fields in which case it is not necessary to solve them explicitly in 
order to get the final form of the constraints. However, it may happen]^ that some 
restrictions on the fields and momenta must be imposed in order to guarantee their 
solvability. These conditions are new secondary constraints that must be put in equal 
footing with the remaining ones, in particular we must impose their stability that 
may lead to new equations and constraints. In the present case a long and tedious 
computation shows that no new constraints must be enforced to ensure the solvability 
of the equations for the Lagrange multipliers. 

By using the explicit expression for the Lagrange multipliers and the fact that, 
once they are plugged in Ht we can express it as the sum of a first class Hamiltonian 
and a linear combination of the primary first class constraints with arbitrary functions 

1. ^- 1 ■ 2 ■ ^ ■ ^ ■ 

of time as coefficients we find out that A*, A^, r/*, 77*, /3* are arbitrary functions 
and the following first class constraints 



2Vair + ^^Jk (k - r'^WVcj) Vak = (27) 



13 



As, for example, if one considers the action (|^). 
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where the last two equations can be substituted by 

2 

Va^f = 

1 

Van = (28) 

by using the first pair. We also find the following second class constraints 

1 

- Bt = 

I • -k = o 

a\ = (29) 



2 



The symplectic structure for the initial set of canonical variables (after the elimination 
of the arbitrary objects A*, A*, r/*, 77*, /?*) is 

n = -2 I d^x [dAl^ix) Adl^^(a;) + dl5,"(s) Adla^a;) + (^hii^) Adl^^ ,''(a;) + 

2 5 2 

+dAi{x) Adll^(x) + dl5,"(x) Adla^x) + dl^i(x) Ad^/i • (x) 
Now we must pull it back onto the hypersurface in phase space defined by the second 

12- ^ 

class constraints. As a first step we use (|29D to remove aa, o\-, -B", and B\. At this 
stage we are left with the first class constraints 

i'^ + lr-r^'^'^'^fe^e. = 

Va^- = (30) 



the second class constraints 



lh-(rr''^m\-^^i)=^. (31) 
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n= -2 d^x 

JT. 



and the symplectic structure 

'^A\{x) Adr^^(a;) + dl^l(x) Adl^/;^(x) + 

+>lA\{x) Adll^(x) + Adl^^(x) 
We have now two possible choices: 

12 1. 2. 

i) Remove and iplhy using (|3TD to get a phase space coordinatized by 
TTa; tTq, "^q, and the first class constraints ( |30D and the symplectic structure 



Vt = - d^x 

JT. 



2^A\{x) Ad."nt{x) + dl^j,(a;) Ad (r^^'^^Vfe^ - |,"(a;) ) + 



1 



+2dl Adllf(x) + dl^i(x) Ad (fi'''"'Wbhci{x) - ^'^[x, 

As we have 54 canonical variables and 24 first class constraints we have three physical 
degrees of freedom per space point. The Dirac brackets of the phase space variables 
can be obtained in closed form by inverting the previous expression for Q in matrix 
form. 

1 2 1.2.1.2_^ 

ii) Remove vr" and vrf to get a phase space coordinatized by A^^, A^^, ipl^, 7]\^ 
and rj\. The first class constraints are given now by 

2 

WaFt = 
1 

"WaFt = 

and the symplectic structure is 

n = -'^J^d^x [dAiix) Ad [fj'^^'WbV ci{x) - 2^'^{x)) + d^Hx) Ad^^(x) + 

+di^(x) Ad {fi'^'^^b^ridx) - 2^^(a;)) + d^^(x) Ad^^(x)" . 

As before we have 54 canonical variables and 24 first class constraints per space point. 
In both cases we can recover in a straightforward way the formulation of ^ by using 
the consistent gauge fixing conditions 
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If instead of using the action @ we take (|^ the procedure is considerably more 
involved due to the appearance of several layers of secondary constraints that come 
up as consistency conditions for the solvability of the Lagrange multipliers equations 
that are avoided here by using the Stiickelberg procedure. Anyway, as argued in the 
main text, the formulation is strictly equivalent to the other ones presented in the 
paper. 
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